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Abstract. The exchange of electroweak gauginos in the f— or m— channel allows squark pair 
production at hadron colliders without color exchange between the squarks. This can give rise to 
events where little or no energy is deposited in the detector between the squark decay products. 
We discuss the potential for detection of such rapidity gap events at the Large Hadron Collider 
(LHC). We present an analysis with full event simulation using PYTHIA as well as Herwig++, 
but without detector simulation. We analyze the transverse energy deposited between the jets from 
squark decay, as well as the probability of finding a third jet in between the two hardest jets. For 
the mSUGRA benchmark point SPSla we find statistically significant evidence for a color singlet 
exchange contribution. 
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SUSY RAPIDITY GAPS 

One of the main objectives of tlie Large Hadron Collider (LHC) is the search for super- 
symmetric (SUSY) particles [1]. In the energy range of the LHC we expect squark pair 
production to be one of the most important channels for the production of superparticles 
[0]. Squark pair production includes contributions with electroweak (EW) exchange par- 
ticles at tree [11,13] or one-loop ^ level. The tree-level EW contributions can change the 
production cross section by up to 50% Moreover, EW gaugino exchange in the t— or 
M— channel gives rise to events with no color connection between the produced squarks. 
QCD radiation then preferentially takes place in the phase space region between the re- 
spective color connected initial quark and final squark, not between the two outgoing 
squarks. If the rapidity region between both squarks is indeed free of QCD radiation 
it is called a "rapidity gap". The situation is different for the lowest order QCD con- 
tribution. The final squarks are color-connected and radiation into the region between 
them is expected. This difference might allow to isolate events with electroweak gaugino 
(color singlet [CS]) exchange, which could e.g. lead to new methods to determine their 
masses and couplings. The above discussion describes a single partonic reaction produc- 
ing stable squarks. In reality, the squarks will decay. Even if we assume that each squark 
decays into a single jet , the rapidity distribution of these jets will differ from that of the 
squarks. Squark decay also leads to additional parton showers from final state radiation. 
Moreover, the underlying event produced by the beam remnants and their interactions 
can also deposit energy in the gap. 
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FIGURE 1. Transverse energy in the rapidity-gap region as predicted by full event simulations of 
squark pair production using HERWIG++ and PYTHIA 6.4. Black histograms show pure QCD results, 
while the red (gray) histograms include electroweak contributions. The errors are statistical only. 



NUMERICAL SIMULATION 

In the following we perform a full event simulation, including squark decays, hadroniza- 
tion, jet reconstruction and the "underlying event", however without simulating the de- 
tector To be specific, we consider rapidity-gap events for squark pair production at the 
LHC, where electroweak (EW) contributions at tree level are included. The production 
of the first two generations of squarks via t— and m— channel diagrams is taken into 
account, j'— channel contributions are neglected, since they are quite small and are 
not expected to lead to rapidity-gap events. The mass spectrum and branching ratios 
of the sparticles are obtained from SPheno [6]. Analytical expressions for the squared 
and averaged matrix elements are given in Ref. [|4t]. We implemented the relevant matrix 
elements for QCD and as well for EW contributions in a simple parton-level simulation. 
Jets are reconstructed via the kj clustering algorithm of FastJet [7]. Events are analyzed 
using the program package root [8]. 

We impose the following cuts: the two highest transverse momentum jets have to 
satisfy: ErUi) > 100 GeV ; | rj {jt) | < 5.0 (z = 1 , 2) . We further suppress SM back- 
grounds by requiring a large amount of missing transverse energy: > 100 GeV. 
Squark pair events containing at least one SU (2) singlet squark, which give rise to quite 
small EW contributions [|4|], are suppressed by requiring the existence of two like-sign 
charged leptons, with pr(^0 > 5 GeV; |t](£,)| < 2.4 (z = 1,2). In order to be able 
to define a meaningful rapidity gap, the two leading jets should be well separated in 
rapidity: At] = \T]{j\) — ^7(72)! > 3.0. We have to take into account that the two jets 
have finite radii. The "gap region" is therefore defined as: min[7](ji),77(j2)] -1-0.7 < 
7] < max [77 (j 1 ) , T] (j2)] — 0.7 . One can expect that most of the particles produced during 
hadronization are within the cone of 0.7 of the corresponding jets [9]. The overall effi- 
ciency of our cuts relative to the entire squark pair sample is less than 1%, cf. [10]. Since 
we want to avoid "event pile-up", i.e. multiple pp interactions during the same bunch 
crossing, we assume an integrated luminosity of 40 fb^^ at ^ = 14 TeV 

Our first attempt to isolate "rapidity gap events" uses a completely inclusive quantity. 
We define ^^f^paj-tides '^he total transverse energy deposited in the gap region; this is 
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computed from all photons and hadrons in the event (after hadronization and decay of 
unstable hadrons), but does not include the leptons produced in and decays. The 
distribution of ^^f^paj-ticies shown in Fig. [T] for Herwig++ and PYTHIA 6.4. In this 
and all following figures, black and red histograms denote pure QCD and QCD+EW 
predictions, respectively. We also show the statistical error for each bin. 

We note that including EW contributions increases the number of events, although 
in most bins this effect is statistically not very significant. However, in the first bin, 
where the total Ej is less than 5 GeV, the inclusion of these CS exchange contributions 
increases the number of events by a factor of 2.8 ± 1.1 and 2.36 ±0.56 in the Herwig++ 
and PYTHIA 6.4 simulations, respectively. This indicates that CS exchange does lead to 
"gap" events where little or no energy is deposited between the two hard jets. 

The difference between the two generators is as large as the effect from the CS 
events: PYTHIA 6.4 without CS exchange contributions predicts almost exactly the 
same number of events in the first bin as Herwig++ with CS exchange. PYTHIA 6.4 
also predicts a ^^f^paj-tides distribution which is quite flat beyond 20 GeV, whereas the 
distribution predicted by Herwig++ flattens out only at about 40 GeV. One might thus be 
able to use the higher bins, where the effect of the CS exchange contributions is not very 
sizable, to decide which generator describes the data better, or to tune the Monte Carlo 
generators to the data. This should reduce the difference between the two predictions. 

Predicting the total transverse energy flow is difficult, since this observable is strongly 
affected by semi- and non-perturbative effects. Thus, we discuss in the following the 
occurrence of relatively soft "minijets" in the gap region. Fig. [2] shows the fraction of 
events where the energy £^f^et most energetic jet in the gap region is less than the 
value £'f^et max displayed on the .x:— axis. Since Fig. |2] shows event fractions, all curves 
asymptotically approach 1 at large £'|'^et ^lax- assume that jets with transverse energy 
above ^r.thresh = 5 GeV can be reconstructed. If the true threshold is higher, the curves 
should simply be replaced by constants for £^f^et max — ^r,thresh- We expect that the 
underlying event by itself generates few, if any, reconstructable jets. The results here 



are nevertheless still not quite immune to non-perturbative effects, since reconstructed 
jets may also contain a few particles stemming from the underlying event. Note also that 
a jet whose axis lies in the gap region might contain (mostly quite soft) particles that 
lie outside of this region. Conversely, even though we use a cluster algorithm, where by 
definition each particle belongs to some jet, some particles in the gap region might be 
assigned to a jet whose axis lies outside this region. 

We see that PYTHIA 6.4 predicts more events without jet in the gap region than 
Herwig++. This is consistent with the observation that PYTHIA predicted more events 
with little or no energy deposited in the gap region and tends to generate fewer hard 
gluons than Herwig++ does, see [10]. Moreover, both PYTHIA 6.4 and Herwig++ 
predict a significant increase of the fraction of events without jet in the gap region once 
EW, CS exchange contributions are included; the effect is statistically most significant 
for £^|^et niax ~ GeV. Here both generators predict an increase of the fraction of 

events without (sufficiently hard) jet in the gap by about 0.05. Taking a threshold energy 
of 30 GeV as an example, Herwig++ predicts about 1,570 out of the total of 4,032 events 
without jet in the gap once EW effects are included. This should allow to measure the 
fraction of events without jet in the gap to an accuracy of about 0.01 ; a shift of the jet-in- 
the gap fraction by 0.05 thus corresponds to a change by about five standard deviations 
(statistical error only). 

Unfortunately a pure SUSY QCD PYTHIA simulation leads to a fraction of events 
without jet of 0.45 for the same threshold energy, which is higher than the Herwig++ 
prediction including EW contributions. Clearly these large discrepancies between the 
two MC generators have to be resolved before reliable conclusions about the color flow 
in squark pair events can be drawn. In Ref. [[lol] we showed that a similar difference 
between the predictions of the two generators also exists for the jet-in-the-gap cross 
section in standard QCD. This indicates that such standard QCD events can be used 
during the earliest phases of LHC running to improve the event generators, hopefully 
to the extent that the difference between their predictions becomes significantly smaller 
than the effect from the EW contributions. 
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